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Abstract 

With regard to homeoidally striated Jacobi ellipsoids (Caimmi &: 
Marmo 2005), a unified theory of systematically rotating and peculiar 
motions is developed, where both real and imaginary rotation (around 
a fixed axis) are considered. The effect of positive or negative resid- 
ual motion excess, is shown to be equivalent to an additional real or 
imaginary rotation, respectively. Then it is realized that a homeoidally 
striated Jacobi ellipsoid with assigned velocity distribution, always ad- 
mits an adjoint configuration i.e. a classical Jacobi ellipsoid of equal 
mass and axes. In addition, further constraints are established on 
the amount of residual velocity anisotropy along the principal axes, 
for triaxial configurations. Special effort is devoted to investigating 
sequences of virial equilibrium configurations in terms of normalized 
parameters, including the effects of both density and velocity profile. 
In particular, it is shown that bifurcation points from axisymmetric 
to triaxial configurations occur as in classical Jacobi ellipsoids, con- 
trary to earlier results (Wiegandt 1982a, b; Caimmi & Marmo 2005). 
The reasons of the above mentioned discrepancy are also explained. 
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The occurrence of centrifugal support at the ends of major equato- 
rial semiaxis, is briefly discussed. An interpretation of recent results 
from numerical simulations on stability (Meza 2002), is provided in 
the light of the model. A physical interpretation of the early Hubble 
sequence is shortly reviewed and discussed from the standpoint of the 
model. More specifically, (i) elliptical galaxies are considered as iso- 
lated systems, and an allowed region within Ellipsoidland (Hunter & 
dc Zccuw 1997), related to the occurrence of bifurcation points from 
ellipsoidal to pear-shaped configurations, is shown to be consistent 
with observations; (ii) elliptical galaxies are considered as embedded 
within dark matter haloes and, under reasonable assumptions, it is 
shown that tidal effects from embedding haloes have little influence 
on the above mentioned results; (iii) dark matter haloes and hosted 
elliptical galaxies, idealized as a single homeoidally striated Jacobi el- 
lipsoid, are considered in connection with the cosmological transition 
from expansion to relaxation, by generalizing an earlier model (Thuan 
Sz Gott 1975), and the existence of a lower limit to the fiattening of 
relaxed (oblate-like) configurations, is established. On the other hand, 
no lower limit is found to the elongation of relaxed (prolate-like) con- 
figurations, and the observed lack of elliptical galaxies more elongated 
than E7 needs a different physical interpretation, such as the occur- 
rence of bending instabilities (Polyachenko & Shukhman 1979; Merritt 
& Hernquist 1991). 

keywords - cosmology: dark matter - galaxies: evolution - galaxies: 
formation - galaxies: haloes - galaxies: structure. 

1 Introduction 

Large-scale celestial objects, represented as self-gravitating fluids, exhibit 
different features according if their subunits, or "particles" , are "coUisional" 
or "collisionless" . 

In the former alternative, the gravitational field that is generated by the 
system as a whole is negligible in respect of the force between two colliding 
subunits, when they repel each other strongly. Consequently, particles in 
self-gravitating, coUisional fluids, are subjected to violent and short-lived 
accelerations as they are sufficiently close to each other, interspersed with 
longer periods when they move at nearly constant velocity. 
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In the latter alternative, the gravitational field that is generated by the 
system as a whole is dominant in respect of the force between two sub- 
units, even if they are close to each other. Consequently, particles in self- 
gravitating, coUisionless fluids, are subjected to smooth and long-lived accel- 
erations through their trajectories. 

Gas in stars and stars in stellar systems may be approximated, to a good 
extent, as coUisional and coUisionless, ideal self-gravitating fluids, respec- 
tively. The statistical problem of particle motion in the latter case is that of 
the former, with the collisions left out. "Ideal" has to be intended as "par- 
ticles collide as perfectly and undeformable spheres" in the former situation, 
and "particles do not interact each other at all" in the latter. 

The motion equation of fluid flow turns out to be the same for both 
coUisional and coUisionless fluids (e.g., Jeans 1929, Chap. II, §26, Chap. VII, 
§§ 211-215; Binney & Tremaine 1987, Chap. 4, § 4.2) provided gases are thought 
of as far from equilibrium. Then the velocity distribution of molecules no 
longer obeys Maxwell distribution and is, in general, anisotropic; accordingly, 
the pressure is represented by a stress tensor. In the special case of isotropic 
velocity distributions, the pressure attains its usual meaning and the motion 
equation reduces to Euler's equation, provided the velocity of an infinitesimal 
fiuid element is intended as the mean velocity of all the particles within the 
same element at the time considered. On the contrary, Euler's equation for 
a given fiuid with isotropic distribution can be generalized to an anisotropic 
one, by replacing velocities with mean velocities and pressures with stress 
tensors, in the sense specified above. 

For special classes of ideal, self-gravitating fiuids, such as steadily ro- 
tating polytropes with polytropic index n > 1/2 (Vandervoort 1980a) and 
isothermal spheres (e.g., Binney & Tremaine 1987, Chap. 4, §4.4b), a one-to- 
one correspondence has been discovered between coUisional and coUisionless 
systems with equal physical parameters. Steadily rotating polytropes may 
be represented, to a first extent, as steadily rotating, homeoidally striated 
ellipsoids (e.g., Vandervoort 1980b; Vandervoort & Welty 1981; Lai et al. 
1993). 

Though most astronomical bodies exhibit ellipsoidal- like shapes and isopy- 
cnic surfaces, still some caution must be used in dealing with the above 
approximation, with regard to local values of physical parameters. This is 
why steadily rotating polytropes are in hydrostatic equilibrium while, ow- 
ing to Hamy's theorem (e.g., Chambat 1994), the contrary holds for their 
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homeoidally striated ellipsoidal counterparts. On the other hand, homcoidally 
striated ellipsoids may safely be assumed as a viable approximation to self- 
gravitating fluids, with regard to typical values of physical parameters, aver- 
aged over the whole volume ^(e.g., Vandervoort 1980b; Vandervoort & Welty 
1981; Lai et al. 1993). 

Celestial objects have been modelled as collisional, self-gravitating fluids, 
in particular homeoidally striated ellipsoids, since the beginning of the sci- 
entific era (see e.g., Chandrasekhar 1969, Chap. 1). The evidence of rotation 
in spiral galaxies and the symmetry of figure shown by elliptical galaxies and 
spiral bulges, suggested the following (e.g., Jeans 1929, Chap. XIII, §299): 
(i) all (regular) galaxies rotate; (ii) the symmetry of figure is precisely such as 
rotation might be expected to produce; (iii) the observed shapes of (regular) 
galaxies can be explained as the figures assumed by masses rotating under 
their own gravitation. In particular, the system of primeval stars must have 
conserved roughly the same shape as the original mass of gas from which it 
emerged (e.g., Blaauw 1965). 

The above outlined ideas were accepted up to about thirthy years ago, 
when observations yielded increasing evidence, that giant elliptical galaxies 
cannot be sustained by systematic rotation (e.g., Bertola & Capaccioli 1975; 
Binney 1976; lUingworth 1977, 1981; Schcchtcr & Gunn 1979). Then giant el- 
liptical galaxies must necessarily be conceived as collisionless, self-gravitating 
fiuids, owing their boundaries to anisotropic peculiar velocity distribution of 
stars (Binney 1976, 1978), with a neghgible contribution from figure rota- 
tion. Accordingly, the shape corresponds - in general - to triaxial ^ellipsoids. 
Owing to high-resolution simulations, the same holds also for (nonbaryonic) 
dark matter haloes hosting galaxies and cluster of galaxies (e.g., Hoeft et al. 
2004; Rasia et al. 2004; Baihn & Steinmetz 2004). 

Isotropic peculiar velocity distributions necessarily imply configurations 
which rotate around the minor axis. Accordingly, empirical evidence of sys- 
tematic rotation around the major axis (e.g., Bertola & Gallctta 1978). in 
absence of tidal potential, makes a signature of the occurrence of anisotropic 
peculiar velocity distribution. 

^An important exception is the energy ratio of rotational to random motions, Ej-ot/Epec 
(Caimmi 1979, 1983). 

^Strictly speaking, all tridimensional bodies may be conceived as triaxial, in particular 
spheres and spheroids. Throughout this paper, "triaxial" shall be intended as denoting 
ellipsoids where the axes are different in length. 
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Though anisotropy in pecuhar velocity distribution is a basic ingredient in 
the description and investigation of stellar systems and hosting dark matter 
haloes, conceived as self-gravitating coUisionless fluids (e.g., Binney 1976, 
1978, 1980), still no attempt (to the knowledge of the author) has been 
made to gain more insight into the intrinsic connection between systematic 
and random motions, and the extent to which they are related in determining 
the figure of equilibrium. The current paper aims to overhelm the above 
mentioned gap, using a theory where systematic and random motions are 
unified. 

To this aim, the effect of anisotropy in peculiar velocity distribution shall 
be related to the effect of systematic (real or imaginary) rotation (Caimmi 
1996b). It will be shown the existence of a one-to-one correspondance be- 
tween configurations characterized by systematic {real) rotation and anisotropic 
peculiar velocity distribution, on one hand, and configurations characterized 
by systematic {real or imaginary; prograde and/or retrograde) rotation and 
isotropic peculiar velocity distribution. 

An alternative explanation for the absence of elliptical galaxies (and non- 
baryonic dark matter haloes) more flattened or elongated than E7 as due 
to bending instabilities, has been suggested long time ago from analytical 
considerations involving homogeneous (oblate and prolate) spheroids (Poly- 
achenko & Shukhman 1979; Fridman & Polyachcnko 1984, Vol.1, Chap. 4, 
Sect. 3.3, see also pp. 313-322; Vol. 2, p. 159) and numerical simulations in- 
volving inhomogeneous (oblate and prolate) spheroids (Merritt & Hernquist 
1991; Merritt & Sellwood 1994). 

In a cosmological scenario (Thuan & Gott 1975), the occurrence of a 
limiting ellipticity in oblate configurations depends on the amount of spin 
growth regardless from the onset of bending instabilities. An interesting 
question could be if a similar conclusion holds for prolate configurations. 
It will be seen that the onset of bending instabilities is necessary for the 
occurrence of a limiting ellipticity in prolate configurations. 

The current paper is organized as follows. The general theory of homcoidally 
striated density profiles, and the results of interest to the aim of this inves- 
tigation, are outlined in Sect. 2. A unified theory of systematic and random 
motions, implying a definition of imaginary rotation, is developed in Sect. 3 
where, in addition, an interpretation is provided to recent results from numer- 
ical simulations on stability. An application of the above mentioned theory, 
namely a physical interpretation of the early Hubble sequence, is performed 
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in Sect. 4, on general grounds, by considering (i) elliptical galaxies as isolated 
systems; (ii) elliptical galaxies as embedded in dark matter haloes; (iii) dark 
matter haloes and hosted eUiptical galaxies, ideahzed as a single homeoidally 
striated Jacobi ellipsoid, in connection with the cosmological transition from 
expansion to relaxation; and establishing the existence of a lower limit to the 
flattening of relaxed (oblatc-likc) configurations. Some concluding remarks 
are drawn in Sect. 5, and a few arguments are treated with more detail in 
the Appendix. 

2 General theory 

2.1 Homeoidally striated density profiles 

A general theory for homeoidally striated density profiles has been developed 
in earlier approaches (Roberts 1962; Caimmi 1993a; Caimmi & Marmo 2003, 
2005 ^, hereafter quoted as CM03, CM05, respectively), and an interested 
reader is addressed therein for deeper insight. What is relevant for the current 
investigation, shall be mentioned and further developed here. 

The isopycnic (i.e. constant density) surfaces are defined by the following 
law: 



where po; ''"o; are a scaling density and a scaling radius, respectively, related 
to a reference isopycnic surface, ^ is a scaled distance, independent of the 
direction along which radial coordinates are calculated, and E — ^{r), is re- 
lated to the boundary, where r = R. Both cored and cuspy density profiles, 
according to the explicit expression chosen for the scaled density, f{C,), are 
represented by Eqs. (1). In the former alternative, it is used a different nor- 
malization with respect to Caimmi (1993a), where ^ = r/R and po is the 
central density. 

The mass, the inertia tensor, and the potential self-energy tensor are: 



more extended version is available at the arxiv electronic site, as astro-pli/505306. 
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Ipq = 6pqUinrMal ] (3) 

\^sel)pq — ^sel\^sel)pq — ^pq i \ ) 

3 Q^2 GM\ 
Esel — Z^\Esel)ii — ^selBsel — S \ (5) 

3 

{Bsei)pq — Spqep2epsAp ; S^e/ = ^ ^82653^5 ; (6) 

Spq = l'sel{Bsel)pq S = VgelBgel ', (7) 

where 6pq is the Kronecker symbol; G is the constant of gravitation; Umas, 
^inr-i ^sei, are profile factors i.e. depend only on the density profile via the 
scaled radius, S; ai, 02, 03, are semiaxes; e^g = ap/aq are axis ratios; Ai, A2, 
As, are shape factors i.e. depend only on the axis ratios; and Mq is the mass 
of a homogeneous ellipsoid with same density and boundary as the reference 
isopycnic surface: 

An 

Mo = — poaoiao2ao3 ; (8) 

where aoi, 002, ^03, are the semiaxes of the ellipsoid bounded by the reference 
isopycnic surface. 

The limiting case of homogeneous configurations reads: 

/(0 = 1 , 0<e<H ; (9a) 

I ,2 X 3 /A'U\ 

^mas — ^ ) ^inr — T j ^sel — TTT > V / 

o iU 

for further details, see CM03, CMOS. 



2.2 Homeoidally striated Jacobi ellipsoids 

In dealing with angular momentum and rotational energy, the preservation 
of (triaxial) ellipsoidal shape imposes severe constraints on the rotational ve- 
locity field. Leaving an exhaustive investigation to more refined approaches, 
our attention shall be hmited here to the special case of homeoidally stri- 
ated Jacobi eUipsoids (CMOS). Accordingly, the systematic velocity field is 
defined by the law: 

Vrot{r;e) ^ Vrotiap,0) _ 

Vrot{R,0) Vrot{ap,0) 
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or equivalently, the angular (with respect to X3 axis) velocity field is defined 
by the law: 

Qjr, e) _ n{a'^. 0) _ 

n{R,e) Q(ap,o) ' ^ ' 

where {r,9), {R,0), represent a point on a generic isopycnic surface and on 
the boundary, respectively, along a fixed radial direction, and (a^, 0), (op, 0), 
p = 1,2, represent the end of the corresponding equatorial semiaxis. 

It is worth noticing that the following special cases are described by 
Eqs. (10) or (11): (a) rigid rotation; (b) constant rotation velocity every- 
where; (c) rigid rotation of isopycnic surfaces and constant rotation velocity 
on the equatorial plane. To maintain ellipsoidal shapes, differential rota- 
tion [e.g., cases (b) and (c) above] must necessarily be restricted to axysim- 
metric figures i.e. spheroidal configurations. In the limiting situation of 
homogeneous, rigidly rotating, dynamical (or hydrostatic) equilibrium con- 
figurations, homeoidally striated Jacobi ellipsoids reduce to classical Jacobi 
ellipsoids (CM05). 

The angular-momentum vector and the rotational-energy tensor are: 

Js^ Ss3r]anmJ^anmMap{l + elp){Vrot)p P Q ^ ; (12) 

{Erot)pq = 5pq0- " ^pi)'nrot^rotM [{Vrotjpf ] (13) 

and the related module and trace, respectively, read: 

el^)ai{vrot)i ; (14) 

Erot = rirotVrotM{l + £21) [{Vrot)lf ] (15) 

where r]anmi^anm, Vrot, are shape factors, z/„„„, Urot, are profile factors, {Vrot)p 
is the rotational velocity at the end of the semiaxis, Op, p = 1,2, and the 
rotation axis has been chosen to be x^. For further details, see CMOS, CMOS. 
The combination of Eqs. (14) and (15) yields: 

iyE^ot)pq — ^ram (-^ram)^? — McL^ '^pq ) (16a) 

; (16b) 



i/2 

anm 



{Bram)pq = 8pq{l - Sps) J° J' ^ ] (I6c) 

nanm \^ ^ ^21) 

T^pq ^rami^ramjpq j (16d) 
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which makes an alternative expression of the rotation-energy tensor, and: 



J2 J2 
Maf Maf 

^ ^ramBram > (■'-'^^) 

which makes an alternative expression of the rotation energy. 
The limiting situations outlined above, read: 

O 

for homogeneous configurations in rigid rotation, according to case (a); 

_ 1 _ 1 _ . 

^anm — ^ j ^rot — g j ^ram — j v / 

for homogeneous configurations with constant velocity on the equatorial 
plane, according to cases (b) and (c); 

■nanm = 1 ; Vrot = ^ ; (20) 

for rigidly rotating isopycnic surfaces, cases (a) and (c); 

Stt 3 . . 

Vanm — j Vrot — ^ j \^^) 

for constant velocity everywhere, case (b). Further details can be found in 
CM03, CMOS. 



2.3 Virial equilibrium configurations 

The virial theorem holds for mass distributions which always remain confined 
within a finite region of space, and the virial equations are satisfied, with 
regard to values of parameters averaged over a sufficiently long time, t ^ T, 
where T is a characteristic period of the system (e.g.. Landau & Lifchitz 
1966, Chap. II, § 10). Let us define virial equihbrium as characterized by the 
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validity of the virial theorem, and relaxed and unrelaxed configurations as 
systems where virial equilibrium does and does not coincide, respectively, 
with dynamical (or hydrostatic) equilibrium (CMOS). 

Let us define systematic and peculiar motions as related to velocity fields 
where the mean value is different from or equal to zero, respectively. For in- 
stance, star motions may safely be conceived as systematic in galactic disks, 
and peculiar in galactic spheroids. In the formulation of the virial equa- 
tions, systematic motions related to axial rotation shall explicitly be treated, 
and the remaining contributions (e.g., vorticity, streaming, peculiar motions) 
shall globally be considered together, and denoted as "residual motions". If 
peculiar motions may safely be thought of as dominant, apart from systematic 
rotation, then in the following "residual" has to be intended as "peculiar" . 

With regard to unrelaxed configurations, the generalized tensor virial 
equations read (CMOS): 

{Eggl^pq + 1{^E.fQt)pq + 2,(^pqEres ~ ] (^2) 

Cp, = ^%^ ; P=l,2,3 ; g = l,2,3 ; (23) 

^res 

ECpp = t^ = C; o<c,,<C; (24) 

p=l ^res 

Cp5 = ; p^q ; (2S) 

where (pq is the generalized anisotropy tensor, E^es is the residual energy, 
and {Eres)pq — CpqEres IS the effective residual-energy tensor i.e. the right 
amount needed for the configuration of interest to be relaxed. The diagonal 
components of the generalized anisotropy tensor, (^pp, may be conceived as 
generalized anisotropy parameters. The related trace, C, may be conceived 
as a virial index, where C = 1 corresponds to null virial excess, 2AEres = 
2{Eres — Eres), which does not necessarily imply a relaxed configuration ^, 
C > 1 to positive virial excess, and C < 1 to negative virial excess. 
Let us define the effective anisotropy tensor, as: 

4. = %^ ; P= 1,2,3 ; g = 1,2,3 ; (26) 



*For instance, a homogeneous sphere undergoing coherent oscillations exhibits C > 1 

at expansion turnover and ^ < 1 at contraction turnover. Then it necessarily exists a 
configuration where C = 1 which, on the other hand, is unrelaxed. 
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ECpp = 1 ; 0<Cpp<l ; (27) 

p=i 

Cp. = ; P^q \ (28) 

where the diagonal components, ^pp, may be conceived as effective anisotropy 
parameters. The combination of Eqs. (23), (24), and (26) yields: 

and the specification of the effective and generalized anisotropy parameters, 
C,pp and Qpp, and the residual kinetic energy, Eres, implies the specification of 
the effective, residual kinetic-energy tensor, {Eres)pq- 

The substitution of Eqs. (4)-(7) and (16)-(17) into the virial equations 
(22), and the particularization to the rotation axis, p — S, allows the following 
expression of the residual energy (CMOS): 

1 GM^ S33 

Eres = 7: -r- ; (30) 

accordingly, the remaining virial equations read (CMOS): 

{C33S,,-C<i<iSss)-2hCs3n,,^0 ; g = l,2 ; (31) 

where h may be conceived as a rotation parameter. It is apparent that 
Eq. (31) admits real solutions provided the following inequalities hold: 

^<|^ ; g = l,2 ; (33) 

<,33 '->33 

which is the natural extension of its counterpart related to axisymmetric, 
relaxed configurations (Wiegandtl982a,b). 

2.4 Rotation and anisotropy parameters 

To get a more evident connection with the centrifugal potential on the bound- 
ary, let us define the rotation parameter (CMOS): 

"-2^ 
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where Q = Q{ai, 0) is the angular velocity at the end of the major equatorial 
semiaxis, denoted as ai, and p is the mean density of the ellipsoid: 

- 3 M 

P = ; (35) 

47r aia2a3 

the above definition of the rotation parameter, v, makes a generalization of 
some special cases mentioned in the literature (e.g., Jeans 1929, Chap. IX, 
§232; Chandrasekhar & Leboviz 1962). 

The combination of Eqs. (14), (32), (34), and (35) yields: 

7 3 2 2 (1 + ^2l)^ /oc\ 

h = T^Vanml'anm^- ; (36) 

^ ^21631 

which links the rotation parameters, h and v. An explicit expression of the 
rotation parameter, h, may directly be obtained from Eq. (31), as: 

h=l^J^^f^^ ; ,= 1,2 ; (37) 

^ {,33 '<^qq 

and the substitution of Eqs. (37) into (36), using (16) and (17), allows the 
explicit expression of the rotation parameter, v, as: 



_ l_^2l£31 1 C33^gg - Cqq^ss _ ^ _ i r, . (oo\ 

•J U + ^21^ Vamn^amn ^33 '<-qq 



which, in the special case of rigidly rotating, homogeneous configurations 
with isotropic peculiar velocity distribution, reduces to a known relation for 
Jacobi ellipsoids and, with the additional demand of axial symmetry, to a 
known relation for MacLaurin spheroids (e.g., Jeans 1929, Chap. VIII, §§189- 
193; Chandrasekhar 1969, Chap. 5, §32, Chap. 6, §39; Caimmi 1996a). 

An explicit expression of anisotropy parameter ratio, can be obtained via 
Eq. (37). The result is (CMOS): 



TZ 

l-2h^ 
^qq 



Cqq_ _ ^qq 

C33 3-^3 

C22 S22 ~ 2/i7?.22 



q=l,2 ; (39) 

(40) 
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and the combination of Eqs. (24) and (39) yields: 

(41) 



C33 ^33 



C s-2hn ' 

which provides an alternative expression to Eqs. (37) and (38), as: 

2 C337^ ^ ' 

_ 1 £21 £31 1 C33»^ ~ C'^33 / 

''-3(l + 6i0^r7Ln^Ln C337^ ' ^ ^ 

that arc equivalent to Eq. (31), and then admit real solutions provided in- 
equality (33) is satisfied. 

Finally, Eqs. (31) may be combined as: 

^11 _ C33'^ll — Cll^33 . (44^ 
'^22 C33'522 — C22'533 

where it can be seen that Eqs. (40) and (44) are changed one into the other, 
by replacing the terms, S^^Cqq/ C,33-i with the terms, 2h7lqq, and vice versa. 
The above results may be reduced to a single statement. 

Theorem 1. Given a homeoidally striated Jacobi ellipsoid, rotating around 
an axis, x^, the relation: 

Xn Sii — Yn 



X22 S22 — Y22 

Xqq — -r^Sss , 2hTZqq ] Q — 1,2 ] 



Yqq — 2hTZqq , y^S^-^ ; Q — 1,2 ; 
(,33 



is symmetric with respect to Xgq and Yqq, being the former tensor related 
to anisotropic residual velocity distribution, and the latter one to systematic 
rotation around the axis, X3, or vice versa. 

In the special case of axisymmetric configurations, the shape factors re- 
lated to equatorial axes coincide, A2 — Ai, as €21 — 1 (e.g., CMOS). Then 
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iSii = 1S22 owing to Eqs. (6), (7), and TZu = 7^22 owing to Eqs. (17), which 
necessarily imply (n = (22, owing to Eq. (44). 

In the general case of triaxial configurations, the contrary holds, A2 7^ ^1, 
as 621 7^ 1, then Su ^ S22 and TZu 7^ 7^.22, then the equality, Cii = C22, via 
Eq. (40), imphes the vahdity of the relation: 

1 «Sii — ^22 , s 

if otherwise, the residual velocity distribution along the equatorial plane ^ is 
anisotropic i.e. (u 7^ C22- The related degeneracy can be removed using an 
additional condition, as it will be shown in the next section. 
The above results may be reduced to a single statement. 

Theorem 2. Isotropic residual velocity distribution along the equatorial 
plane, Cii = C22; niakes a necessary condition for homeoidally striated Ja- 
cobi ellipsoids to be symmetric with respect to the rotation axis, x^, i.e. 
reduce to homeoidally striated MacLaurin spheroids. 

3 A unified theory of systematic and random 
motions 

3.1 Imaginary rotation 

A unified theory of systematic and random motions is allowed, taking into 
consideration imaginary rotation. It has been shown above that Eq. (31), or 
equivalently one among (37), (38), (42), (43), admits real solutions provided 
inequality (33) is satisfied. If otherwise, the rotation parameter - let it be h 
or V - has necessarily to be negative, which implies, via Eqs. (32) or (34), an 
imaginary angular velocity, iQ, where i is the imaginary unit. Accordingly, 
the centrifugal potential takes the general expression: 

T^{xi,X2,X3) = T^[Q{xi,X2,X3)]'^w'^ ; w'^ = xl + xl ] (46) 

■''Throughout this paper, "along the equatorial plane" has to be intended as "along any 
direction parallel to the equatorial plane" . 
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where the minus and the plus correspond to imaginary and real rotation, 
respectively. The centrifugal force related to real rotation, dT^ /dw, has 
opposite sign with respect to the gravitational force, dV/dw. On the other 
hand, the centrifugal force related to imaginary rotation, dT~ /dw, has equal 
sign with respect to the gravitational force, dV/dw. Then the net effect 
of real rotation is flattening, while the net effect of imaginary rotation is 
elongation, with respect to the rotation axis (Caimmi 1996b). 

To get further insight, let us particularize Eq. (31) to the special case of 
null rotation (h — O). The result is: 

^ = 1^ ; h = ■ q=l,2 ■ (47) 

where the right-hand side, via Eqs. (6) and (7), depends on the axis ratios 
only. It can be seen (Appendix A) that Sqq/Sss > 1 for oblate-like con- 
figurations (ai > 02 > 03), and Sgq/S^s < 1 for prolate- like configurations 
(a2 < Oi < 03), which implies (gq/Css ^ 1 for oblate- like configurations, 
and C99/C33 < 1 for prolate-hke configurations. Accordingly, the net effect 
of positive or negative residual motion excess along the equatorial plane is 
flattening or elongation, respectively. In what follows, it shall be intended 
that residual motion excess is related to the equatorial plane. 

The above considerations show a strict connection between real rotation 
and positive residual motion excess, on one hand, and between imaginary 
rotation and negative residual motion excess, on the other hand. At this 
point, it is natural to ask if (positive or negative) residual motion excess can 
be interpreted as (real or imaginary) rotation, and vice versa. 

3.2 Residual motion excess and rotation 

With regard to the rotation parameter, v, it is convenient to use the more 
compact notation ^: 

vn = V ; (48) 

and the substitution of Eqs. (6), (7), (16), (17), (29), into (38) and (43), 
yields the equivalent expressions: 



VN^Aq-^ejqAs ? = 1, 2 ; (49) 

^,33 



^The factor, 3, does not appear in CM05: vn = 3(fAr)cM05- 
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33 



(50) 



which shows that, for homeoidally striated Jacobi elhpsoids, the normahzed 
rotation parameter, vn, depends on the axis ratios, £21, £31, and a single 
anisotropy parameter, (^33 = C33/C- the special case of homogeneous, 
rigidly rotating configurations, owing to Eqs. (9b), (18), and (20), Eq. (48) 
reduces to: = v. 

In the limit of isotropic residual velocity distribution, Cn — C22 — C33 — 
C/3, Eqs. (49) and (50) take the simpler form: 

{vN)iso ^ Ag- ejgAs ; 5 = 1, 2 ; (51) 

{vN)iso = Ui + el^A2 - 2el^A3\ ; (52) 

i + £21 L 

where the index, iso, means isotropic residual velocity distribution. 
Accordingly, Eqs. (49) and (50) may be expressed as: 

Vn = {vN)iso - {vN)ani ] (53) 
{vN)am = - l) 4^3 i (54) 

{vN)am = f - 3) ; (55) 

where {vi\i)ani > for oblate-likc configurations, C^q/Css ^ Ij {tJN)ani ^ 
for prolate-like configurations, Cqq/Css ^ 1; ^^'^ the index, ani, means con- 
tribution from residual motion excess. Accordingly, positive residual motion 
excess along the equatorial plane is related to real rotation. On the con- 
trary, negative residual motion excess along the rotation axis is related to 
imaginary rotation. 

Let us rewrite Eq. (53) as: 

ivN)iso ^Vn+ {vN)ani ] (56) 

which, owing to Eqs. (34) and (48), is equivalent to: 

^2 ^02 ^2. . (57) 
ISO ' am ' v"'/ 
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where the plus corresponds to real rotation, ('UAr)am > 0, and the minus to 
imaginary rotation, {vN)am < 0. Then the effect of residual motion excess 
on the shape of the system, is virtually indistinguishible from the effect of 
additional systematic rotation, with similar (rotation) velocity distribution as 
the pre-existing one. The related, explicit expression, according to Eqs. (11), 
(34), and (55), is: 

^Ini = ^Iniiai) = T{vN)ani2nGp ; (59) 

where the double sign ensures a non negative value on the right hand-side 
member of Eq. (59). The above results may be reduced to a single statement. 

Theorem 3. Given a homeoidally striated Jacobi ellipsoid, the effect of 
(positive or negative) residual motion excess is equivalent to an additional 
(real or imaginary) rotation, with equal spatial distribution as the pre-existing 
rotation, but value at the end of major equatorial axis conform to the rela- 
tion: 

T^lm = Tf^L(«i) = 27rGp f f - 3) ; (60) 

with regard to an adjoint configuration where the residual velocity distribu- 
tion is isotropic. 

Accordingly, a homeoidally striated Jacobi ellipsoid with assigned system- 
atic rotation and residual velocity distribution, with regard to the shape, is 
virtually indistinguishable from an adjoint configuration of equal mass and 
axes, systematic rotation velocity distribution deduced from Eqs. (11), (58) 
and (59), and isotropic residual velocity distribution. 

Owing to Eqs. (53), (54), and (55), Eqs. (51) and (52) are equivalent to 
(49) and (50), respectively. On the other hand, Eqs. (51) and (52) are vahd 
for classical Jacobi ellipsoids with rotation parameter, v — {vN)iso- The 
above results may be reduced to a single statement. 

Theorem 4. Given a homeoidally striated Jacobi ellipsoid, the adjoint con- 
figuration, characterized by isotropic residual velocity distribution and nor- 
malized rotation parameter, ('UAr)iso> exhibits same shape as a classical Jacobi 
ellipsoid of equal mass and axes, and rotation parameter, v = {vN)iso- 
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Accordingly, a homeoidally striated Jacobi ellipsoid (in general, an unre- 
laxed configuration) with assigned systematic rotation and residual velocity 
distribution, with regard to the shape, is virtually indistinguishable from an 
adjoint, classical Jacobi ellipsoid (a relaxed configuration) of equal mass and 
axes, and rotation parameter equal to the normalized rotation parameter of 
the original configuration. 



3.3 Axis ratios and anisotropy parameters 

The combination of alternative expressions of the rotation parameter, {vN)iso, 
expressed by Eqs. (51), yields: 

el,{A2-A^)^{l-el,)el,As ; (61) 

which, for axisymmetric configurations, reduces to an indeterminate form, 
= 0. 

The combination of alternative expressions of the rotation parameter, 
{yN)ani, expressed by Eqs. (54), yields: 

C33-C22 = e^i(C33-Cii) ; (62) 

which, for isotropic residual velocity distributions, reduces to an indetermi- 
nate form, = 0. In addition, axisymmetric configurations (e2i = 1) nec- 
essarily imply isotropic residual velocity distributions along the equatorial 
plane, Cii = (22- 

The combination of Eqs. (24) and (62) yields: 



C - (2 - elX: 



Cii = ^^=-^^ ; (63a) 
^" ^21 



4iC + (1 - 2el)C: 



C22 - ' ;i r''^'' ; (63b) 
J- + ^21 

which, for axisymmetric configurations (e2i = 1) reduces to Cii = C22 = 
(C - (33)72, and the special case, C33 = C/3, reads Cii = C22 = C/3- 

The conditions, Cii > 0, related to Eqs. (23) and (24), and < 621 < 1, 
necessarily imply < C33/C < 1/2- If otherwise, sequences of virial equihb- 
rium configurations cannot attain the oblong shape, 621 — €31 = 0, but must 
stop earlier, where the equivalent relations: 

2 2^33 — C An fRA\ 

£21 = — 7 ; Cii = ; (64) 

S,33 
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are satisfied. The above results may be reduced to the following statements. 

Theorem 5. Given a homeoidally striated Jacobi ellipsoid, the anisotropy 
parameters along the equatorial plane depend on the equatorial axis ratio 
for triaxial configurations. The related expressions coincide, (^u = (,22, In the 
limit of axisymmetric configurations, £21 = 1. 

Theorem 6. A necessary and sufEcient condition for a homeoidally striated 
Jacobi ellipsoid to have isotropic residual velocity distribution, is that the 
anisotropy parameter along the rotation axis, equals one third. 

Theorem 7. Given a sequence of homeoidally striated Jacobi ellipsoids, the 
ending point occurs when the equatorial axis ratio is null, 631 = 0, and/or the 
generalized anisotropy parameter related to the major axis is null, (n = 0, 
which is equivalent to €21 = (2 — C/Css)^^^- The special case of dynamical (or 
hydrostatic) equilibrium, C — implies centrifugal support along the major 
axis provided Cii = 0. 

Accordingly, with regard to homeoidally striated Jacobi ellipsoids, the 
anisotropy parameters along the equatorial plane, Cii and (22, cannot be 
arbitrarily assigned, but depend on the equatorial axis ratio, £21, conform to 
Eqs. (63). On the other hand, the further knowledge of the meridional axis 
ratio, 631, and the rotation parameter, vj\f, allows the determination of the 
rotation parameter, {vj^)ani, via Eqs. (51), (52), (53), and then the ratios, 
Cgg/CsS) C33/C) via Eqs. (54), (55), respectively, or the anisotropy parameter 
along the rotation axis, C33, provided the virial index, defined by Eq. (24), 
is assigned. 

In conclusion, with regard to homeoidally striated Jacobi ellipsoids de- 
fined by assigned axis ratios, eai, €21, rotation parameter, vn, and virial 
index, (, the anisotropy parameters, Cii, C22, C33) cannot arbitrarily be fixed, 
but must be determined as shown above. 

3.4 Sequences of virial equilibrium configurations 

With regard to homeoidally striated Jacobi ellipsoids, it has been shown 
above that adjoints configurations are characterized by (i) centrifugal po- 
tential, Tiso{xi,X2,X3) = T{xi, X2, X3) + Tani{xi, X2, X3), Or Qf^^{xi,X2,Xs) = 
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Q'^{xi,X2,Xs) =F ^araj(^i5 ^2, a^s), where the angular velocity distribution re- 
lated to residual motion excess, is defined by Eqs. (58) and (60); and (ii) 
isotropic residual velocity distribution. Owing to Theorem 3, a sequence of 
homeoidally striated Jacobi ellipsoids in virial equilibrium coincides with the 
sequence of adjoints configurations. The last, in turn, owing to Theorem 
4, coincides with the sequence of classical Jacobi ellipsoids provided the re- 
lated rotation parameters coincide, v = (i^7v)iso; conform to Eqs. (51), (52), 
and imaginary rotation is also considered. Given a homeoidally striated Ja- 
cobi ellipsoid with fixed axis ratios, £21, £31, rotation parameter, {vN)ani, and 
virial index, (, the anisotropy parameters, (u, C22, Cssi are determined via 
Eqs. (51)-(55) and (63). Negative values of the rotation parameter, {vN)iso, 
extend the sequence of classical MacLaurin spheroids to imaginary rotation 
i.e. prolate configurations where the major axis coincides with the rotation 
axis. 

Once more owing to Theorem 4, the bifurcation point of a sequence of 
homeoidally striated Jacobi ellipsoids in virial equilibrium coincides with its 
counterpart along the sequence of classical Jacobi ellipsoids, [631, {vN)iso] = 
(0.582724, 0.187115) (e.g., Jeans 1929, Chap. VIII, § 192; Chandrasekhar 1969, 
Chap. 5, § 33). Accordingly, the axis ratio of the configuration where the bi- 
furcation occurs, is a solution of the transcendental equation (Caimmi 1996a): 



A, ~5-4eii 
where, owing to Eqs. (118): 



621 = 1 ; (65) 



^<1 ; e3i<l ; (66a) 



41^3 



>1 ; e3i>l ; (66b) 



and further analysis of Eq. (65) shows that two solutions exist. One is related 
to spherical configurations (esi = £21 = 1), which belong to the axisymmetric 
sequence (oi = 02 7^ 03); for this reason, it makes a parasite solution to the 
problem under discussion. On the contrary, the remaining solution defines 
the axisymmetric configuration where bifurcation takes place. 

Sequences of homeoidally striated Jacobi ellipsoids can be deduced from 
Figs. 1 and 2. The sequence of classical MacLaurin spheroids (including 
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Jacobi ellipsoids), extended to the case of imaginary rotation, is shown in 
Fig. 1. The rotation parameter, (i^7v)am, is plotted in Fig. 2 as a function of 
the axis ratio, £31, for different values of the effective anisotropy parameter, 
C33 — C33/C) labelled on each curve. With regard to a fixed anisotropy 
parameter, ^33, the related sequence starts from a non rotating configura- 
tion, Vn = i.e. {vi^)ani = {vN)iso, ^nd cuds at a configuration where 
either £31 = or = 0. The locus of the ending points related to the 
latter alternative, is represented as a long-dashed curve. The sequence can- 
not be continued on the right, as imaginary rotation i.e. larger ^33 would be 
needed and a different sequence should be considered. The effect of positive 
(C33 < 1/3) or negative ((^33 > 1/3) residual motion excess is equivalent to 
an additional real or imaginary rotation, respectively. The horizontal non 
negative semiaxis, £31 > 0, (t'Ar)am — 0, is the locus of configurations with 
isotropic residual velocity distribution, ^33 = 1/3. The vertical straight line, 
£31 = 1, is the locus of spherical configurations. 

With regard to real rotation, the ending point of a sequence occurs at 
either null meridional axis ratio, £31 = 0, and/or null effective anisotropy pa- 
rameter related to major equatorial axis direction, Cn = 0. The occurrence 
of the latter situation is marked by a square on the corresponding sequence 
plotted in Fig. 2, and the related locus is represented as a long-dashed curve. 
Configurations on the left are forbidden, as they would imply negative resid- 
ual energy tensor component, {Eres)ii < 0, to satisfy the virial equation 
(22), which would imply imaginary rotation around major equatorial axis 
i.e. systematic motions other than rotation around the minor axis. 

With regard to imaginary rotation, no ending point occurs and the system 
is allowed to attain negative infinite rotation parameter, {vN)iso — ^ —00, and 
infinite meridional axis ratio, £31 — > +00. The related configuration is either 
a spindle (oi — a2 — 0) or an infinitely high cylinder (ai = 02 < 03 ^ +00) . 

Further inspection to Fig. 2 shows additional features, namely: (i) null left 
first derivatives on each sequence at bifurcation point (not marked therein 
for sake of clarity), and (ii) occurrence of symmetric sequences with respect 
to the horizontal axis (including also forbidden configurations). 

The existence of null left first derivatives on each sequence, at bifurcation 
point, cannot be deduced analytically but, on the other hand, it may be in- 
ferred both numerically and graphycally. Further support to this idea comes 
from the following considerations. Sequences branching from the bifurcation 
point on the left, such as those plotted in Fig. 2, are related to triaxial config- 
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Figure 1: The sequence of classical MacLaurin spheroids, from the start- 
ing point (square) to the bifurcation point (triangle), and classical Jacobi 
ellipsoids, from the starting point (triangle) to the bifurcation point (lower 
latin cross), extended to the case of imaginary rotation (negative values of 
the rotation parameter, v, which equals the normalized rotation parameter, 
{v]^)iso-, related to homeoidally striated Jacobi ellipsoids). Both sequences are 
continued in the region of instability. With regard to MacLaurin spheroids, 
it is also marked: the extremum point of maximum (diamond); the point at 
the onset of dynamical instability (St. Andreas cross); and the point at the 
onset of instability towards pear-shaped configurations (upper latin cross), 
which has its counterpart on the Jacobi sequence. Numerical values are taken 
from Chandrasekhar (1969, Chap. 5, §§32-33; Chap. 6, §§39-41). The hori- 
zontal line, {vN)iso = 0, €31 > 0, is tm? locus of non rotating configurations. 
The vertical hne, £31 = 1, is the locus of spherical configurations. The above 
mentioned lines divide the non negative semiplane, €31 > 0, into four regions 
(from top left in clockwise sense): A - oblate-like shapes with real rotation; 
B - prolate-like shapes with real rotation; C - prolate-like shapes with imag- 
inary rotation; D - oblate-like shapes with imaginary rotation. Regions B 
and D are forbidden to sequences of homeoidally striated Jacobi eUipsoids. 



0.0 



0.5 1.0 
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Figure 2: The normalized rotation parameter, {vN)ani, related to residual 
motion excess, as a function of the axis ratio, £31, with regard to homeoidally 
striated Jacobi eUipsoids. Each curve is labelled by the value of the gen- 
eralized anisotropy parameter, ^33 = C33/C- The horizontal non negative 
semiaxis, €31 > 0, {vN)ani = 0, is the locus of configurations with isotropic 
residual velocity distribution, (^33 = 1/3. The vertical straight line, 631 = 1, 
is the locus of spherical configurations. The generic sequence starts from a 
non rotating configuration (short-dashed curves) and ends at a configuration 
where either £31 = or C33 = (long-dashed curve). The regions above 
the upper short-dashed curve and below the long-dashed curve, respectively, 
are forbidden. The non negative vertical semiaxis, (t^7v)am > 0, 631 = 0, 
corresponds to flat ((^33 = 0) or oblong ((^33 = (22 = 0) conflgurations. The 
curves arc symmetric with respect t<?xhe horizontal axis, until the limiting 
curve, C33 = C = 1) is reached in region (D). The limiting configuration where 
Cii = 0, is marked by a square: no configuration in virial equilibrium is al- 
lowed on the left, as it would imply (u < 0. No configuration is allowed on 
the right of the starting point, as it would imply < 0. The sequence start- 
ing from the bifurcation point (not shown) corresponds to C33 = 0.241648, 
Cii = C22 = 0.379176. 



urations where the equatorial axis ratio does not exceed unity, £21 < 1. The 
above mentioned sequence can also be extended to adjoint configurations 
where €21 — > 1, €21 > 1, which corresponds to replacing each equatorial axis 
with the other. Accordingly, the extended sequence of triaxial configurations 
is symmetric with respect to the axis, 631 = (631)5^^ i.e. the value related to 
the bifurcation point, which implies the occurrence of an extremum point of 
maximum or minimum, and then a null first derivative. The above results 
may be reduced to a single statement. 

Theorem 8. Given a sequence of honicoidally striated Jacobi ellipsoids, the 
contribution of residual motion excess, {vN)ani, to the rotation parameter, 
{'<JN)iso, has a null left first derivative at the bifurcation point, 
[d(i;jv)am/de3i](,3^)-^ = 0. 

The occurrence of symmetric sequences with respect to the horizontal 
axis, is deduced from Eq. (55) using the condition: 

(a3)-'-3 = -(C33)-' + 3 ; (67) 

where C33 = C33/C is related to curves characterized by negative (C33 = C33 ^ 
1/3) and positive (^33 = C33 < 1/3) values, respectively, of the rotation 
parameter, {v]\f)ani, see Fig. 2. Couples of symmetric sequences start from 
(CsaiC^s) = (V3il/3); where each curve coincides with the other, and end 
at (1,1/5). Sequences in the range, < C33 < 1/5, have no symmetric 
counterpart. 

Let a point, P[e3i, {vN)iso], be fixed on the extended sequence of classical 
Jacobi ellipsoids, and the axis ratio, 621, be determined by use of Eqs. (51), 
(52). Let a point, P'le-n, {vN)ani], be fixed on the_ plane, [Oe3i{vN)ani], and 
the generalized anisotropy parameters, (n, (22, Css, be determined by use 
of Eqs. (54), (55), (63). Finally, let the normalized rotation parameter, v^, 
be determined by use of Eq. (53). FoHowing the above procedure, sequences 
of homeoidally striated Jacobi ellipsoids may be generated. For assigned 
density profiles and systematic rotation velocity fields, there are three inde- 
pendent parameters, which may be chosen as two axis ratios, €21, €31, and 
one generalized anisotropy parameter, (^33. 
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3.5 Centrifugal support along the major equatorial axis 

The calculation of the gravitational force, induced by homeoidally striated 
Jacobi ellipsoids, involves numerical integrations (e.g., Chandrasekhar 1969, 
Chap. 3, § 20) and is outside the aim of the current paper. With regard to an 
assigned isopyenic surface internal to, or coinciding with, the boundary, the 
gravitational force, Fq, induced at the end of the related major equatorial 
semiaxis, a'l, satisfies the inequality: 

[FG{a[, 0, 0)]sph < FG{a[, 0, 0) < [Faia,, 0, O)]foc ; 

where the left-hand side is related to spherical isopycnics with unchanged 
major equatorial axes, and the right-hand side to confocal isopyenic surfaces 
from the one under consideration to the centre. Owing to Newton's and 
MacLaurin's theorem (e.g., Caimmi 2003), the following relations hold: 

FG{a[, 0, 0) = -27rGpMMM ; (68) 

where p{a[) is the mean density within the isopyenic surface, and w — 
sph, foe, for the striated sphere and the focaloidally striated ellipsoid sur- 
rounded by the homeoidally striated corona, respectively. 

The balance between gravitational and centrifugal force at the end of the 
major equatorial axis of the isopyenic surface under discussion, reads: 

-27!-Gpy,{a[){Ai),uia[ + Q^^^{a[)a[ = ; (69) 

where the angular velocity, Qiso, takes into account both systematic rotation 
and residual motion excess, according to Eq. (57). On the other hand, the 
generalization of Eq. (34) to a generic isopyenic surface within the boundary, 
reads: 

- ; (70) 

and the combination of Eqs. (48), (56), (57), (69), (70), yields: 

{[viso{a[)UU - ^^l#TT = (AiU ; (71) 

where the index, eq, denotes the rotation parameter, Visoia'i), related to 
centrifugal support at the end of major equatorial axis of the isopyenic surface 
under consideration. The above results allow the validity of the relation: 

{[Viso{a'l)]eq}cof < [Viso{a[)]eq < {[Viso{a[)]eq} sph ] (72) 
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where [viso{ai)]eq is the critical value related to the homeoidally striated 
Jacobi ellipsoid, with regard to an assigned isopycnic surface internal to, 
or coinciding with, the boundary. Then Viso{a[) > {[viso{0'i)]eq}cof and 
i^isoid'i) ^ {[i^iso{(ii)]eg}sph make a sufficient and a necessary condition, re- 
spectively, for the occurrence of centrifugal support at the end of major 
equatorial axis in the case Tinder consideration. 

Owing to Newton, s and MacLaurin's theorem, the condition of centrifugal 
support, Eq. (71), for focaloidally striated ellipsoids surrounded by homeoidally 
striated coronae, coincides with its counterpart related to homogeneous el- 
lipsoids with equal mean density, axis ratios, and velocity field. On the other 
hand, the rotation parameter, Viso{a'i), in the special case of focaloidally stri- 
ated ellipsoids, is also expressed by Eq. (71), particularized to classical Jacobi 
ellipsoids. A comparison with the last part of Eq. (71) shows that centrifugal 
support in focaloidally striated ellipsoids occurs only for flat configurations, 
£31 = 0. Accordingly, Eq. (72) reduces to: 



owing to Eq. (71) and {Ai)sph = 2/3 (e.g., CM03). 

With regard to rigid rotation, it is a well known result that the occurrence 
of centrifugal support, at the end of the major equatorial axis, depends on 
the steepness of the density profile (e.g., Jeans 1929, Chap. IX, §§230-240). 
More precisely, a steeper density profile implies an earlier centrifugal support 
and vice versa, unless the bifurcation point from ellipdoidal to pear-shaped 
configurations is attained, which occurs for nearly homogeneous matter dis- 
tributions. A similar trend is expected to hold for any velocity profile of the 
kind considered in the current paper. 

The combination of Eq. (72) with its counterpart related to the boundary, 
yields: 




(73) 




(74) 



a further restriction to mass distributions obeyng the following law: 




< A; < 3 



(75) 
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where k = 3,1,0, represent homogeneous, isothermal, and Roche-hke mass 
distributions, make Eq. (74) reduce to: 

\aij n^{ai) 

where the last factor equals unity for rigid rotation and the ratio, {ai/a[y, for 
constant velocity along the major equatorial axis. Accordingly, centrifugal 
support at the end of major equatorial axis is first attained on the boundary 
in the former alternative, Viso{o,i) < Viso{(ii)- 
In the latter alternative, Eq. (76) reduces to: 

Visoia'i) = Visa (^^^ ; (77) 

which, for sufficiently steep density profiles, < k < 1, shows a similar trend 
with respect to the former alternative. On the other hand, centrifugal support 
is first attained everywhere along the major equatorial axis for isothermal 
mass distributions, k = 1, and at the centre for sufficiently mild density 
profiles, 1 < A; < 3, which implies Viso{a[) > Viso{(ii). 
The combination of Eqs. (48) and (71) yields: 

{[ivN)isoia[)UU^^^^^^{A,)^ ; (78) 



and Eq. (73) reads: 



< [(t;^).4«'i)]e, < ^^^^^ ; (79) 



in terms of the normalized rotation parameter, v^. 



3.6 Bifurcation points 

Owing to Theorem 4, a bifurcation point from axisymmetric to triaxial con- 
figurations is related to an axis ratio, 631, which is independent of the amount 
of systematic rotation and residual motion excess, according to Eq. (56). To 
gain more insight, let us equalize the alternative expressions of Eq. (49). The 
result is: 

^1 - 7^4^ = ^2 - ^e^2^ ; (80) 

S33 S,33 
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and the combination of Eqs. (63) and (80) yields: 



= +^3 ; (81) 
£21 £21 



and the bifurcation points occur at a configuration, where the axis ratio, 631, 
is the solution of the transcendental equation (Caimmi 1996a): 



lini 44^^ = e^,A3 ; (82) 

£21— »1 i — £21 



where Eqs. (81) and (82) coincide with their counterparts related to isotropic 
residual velocity distributions (Caimmi 1996a,b). The above results may be 
reduced to a single statement. 

Theorem 9. Given a sequence of homeoidally striated Jacohi ellipsoids, the 
bifurcation point occurs at the same configuration, as in the sequence of 
classical Jacobi ellipsoids. 

The contradiction with earlier results (Wiegandt 1982a,b; CMOS) is explained 
in the following way. Let us suppose that the generalized anisotropy parame- 
ters, Cii, C22, CsS) can be arbitrarily fixed regardless from Eqs. (63), and assume 
Cii — C22 also for triaxial configurations. Accordingly, Eq. (80) reads: 

J- — f 21 <i,33 ^21 

and the bifurcation points occur at a configuration where the axis ratio, £31, 
is the solution of the transcendental equation: 

^211 — = i^e^i^s ; (84) 

£21^1 1 — £91 C.'^.'? 



-21 S33 



which coincides with Wiegandt (1982b) criterion for bifurcation, with re- 
gard to homeoidally striated Jacobi ellipsoids in rigid rotation. For a formal 
demonstration, sec Appendix B. Then Wicgandt's criterion for bifurcation, 
expressed by Eq. (84), is in contradiction with Eqs. (63), contrary to the cur- 
rent one, expressed by Eq. (82). 
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3.7 Interpretation of results from numerical simula- 
tions on stability 

Given a homeoidally striated Jacobi ellipsoid, the occurrence of anisotropic 
residual velocity distribution may be related to residual motion excess, via 
Eqs. (54) and (55). In fact, the changes: 



do convert (positive or negative) peculiar motion excess into systematic (real 
or imaginary) rms rotation, and the diagonal components of the residual- 
energy tensor read (^res)pp ^ (£'res)pp - i^Eres)pp = (^res)33- 

The velocity field related to rms rotation has to coincide with the pre- 
existing one, apart from a renormalization, defined by Eq. (59), and the 
mean rotational velocity, Qanii^i, X2,X3){xi + 2:2)''^^^, has necessarily to be 
null within any infinitesimal mass element. 

In this view, recent results from numerical simulations on the stability 
of rotating stellar systems (Meza 2002), find a natural interpretation. Each 
computer run therein proceeds along the following steps. 

(i) Take a spherical system with isotropic peculiar velocity distribution and 

null systematic rotation velocity field. 

(ii) Fix a rotation axis, let it be x^. 

(iii) Reverse a selected sense of rotation of a given fraction of the parti- 
cles, by changing the tangential (parallel to equatorial plane) velocity 
component, v^, into its opposite, —V(j). 

(iv) Make the system evolve from the above defined initial configuration. 

It is worth noticing (Meza 2002) that the distribution function is independent 
of the sign of v^, and the whole set of possible initial configurations under 
discussion, are characterized by the same total kinetic and potential energy. 

The simulations show that spherical systems of the kind considered are 
dynamically stable, irrispective of their degree of rotation. Even if very 
rapidly rotating, boundaries fail to become oblate-like. 





(85c) 
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In the light of the current theory, reversing an assigned sense of rotation of 
a given fraction of the particles, implies a conversion of a fraction of residual 
energy into rotational energy, as: 

/^Eres = -^Erot = -nm < > ; (86) 

<<>-ltW]' ; 0<n<f ; (87) 



N-2n AT 

E ^f-^ ; ; (88) 

1=1 

where n is the number of particles with reversed tangential (parallel to the 
equatorial plane) velocity component, is the total number of particles, and 
m is the mass, assumed to be equal for each particle, as done in numerical 
simulations (Meza 2002). The validity of Eq. (88) does not necessarily rule 
out streaming motions other than systematic rotation. On the contrary, 
streaming motions are allowed provided the existence of a given flow implies 
a counter-flow, such that Eq. (88) is satisfied. 

Before reverting the sense of rotation, the following relations hold in con- 
sequence of isotropic residual velocity distribution: 

< vl >=< vl >=< vj > ; (89a) 
<vl> + <vl>=<vf> + <vl> ; (89b) 

where Vp, p — 1,2,3, and Vy^, w — {xf + x^Y^^, are velocity components 
parallel to coordinate axes and equatorial plane, respectively. After reversion 
of the sense of rotation, the following changes are performed: 

Eres Kes - Um < vl > ] (90) 

- nm <vl> ; (91) 

{Eres)qq ^ {Eres)qq " < > ; Q = 1,2 ] (92) 

Erot -^0 + nm<vl> ; (93) 
iErot)qq^O+l-nm<vl> ] ? = 1, 2 ; (94) 



while the contributions from residual motions along the equatorial plane, 
{Eres)w, and the rotation axis, {Eres)33, remain unchanged. 
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The system being relaxed, = 1, in the case under discussion, the gen- 
erahzed and effective anisotropy parameters, Cpp and Cpp, coincide with their 
usual counterparts, and Eqs. (23) and (26) take the explicit form: 

[ryres > — nm < v^> 

- (E.e.) -nm<vl> ' 

where the special case, n = 0, turns back to isotropic residual velocity dis- 
tributions, Cii = C22 = C33 = 1/3- 

In the extreme case where the sense of rotation is completely reversed, 
n — N/2, the changes expressed by Eqs. (90)- (94) reduce to: 

N 2 

Eres Eres " y"^ < >= -Eres (96) 

{Eres)<t> {Eres)<j> - < t>J >= ; (97) 

{Eres)gq ^ {Eres)gg " < "^J >= ^^res ; g = 1, 2 ; (98) 

N ^ 1 

Erot ^ + —m < v; >= -Eres (99) 

{Erot)gg^O + jm<vl>^^Eres ? = 1, 2 ; (100) 

and the anisotropy parameters, expressed by Eqs. (95), reduce to: 

(i/6)£;,e. 1 ^ ^ . 

_(im^_i . .^Q^,. 

which describes the maximum amount of residual motion excess, related to 
complete reversion of the sense of the rotation. 

Turning back to the general case, a comparison between the expression 
of rotational energy, Eq. (15), with the special one related to the case under 
discussion, Eq. (93), reads: 

rirot^rotM{l + e^i) < vl >a,= UTTl < vl > ] (102) 
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where < f ^ >a^ is the rms systematic rotation velocity at the end of the major 
equatorial 

axis , w — {x\ + x^^/'^ — Oi- The related rms angular velocity is: 

2 

(103) 

of 

and the combination of Eqs. (5), (6), (22), (34), (35), (48), (88), (102), and 
(103) yields after some algebra: 

where E — Egeif^ is the total energy, according to the virial theorem. On 
the other hand, the combination of Eqs. (55), (88), and (95) yields: 



3nm <vl> el.As 4^AE^ _ 



where E = —Eres is the total energy, according to the virial theorem. 

The validity of the above results is related to the assumption of relaxed 
configurations, which implies a virial index, ^ = 1, according to Eqs. (22), 
(23), and (24). In addition, the assumption of sphericity necessarily restricts 
the validity of Eqs. (104) and (105) to spherical systems, where €21 = 631 = 1, 
Ai = A2 = As = 2/3, Bsei = 2. Accordingly, the expression of the rotation 
parameters, vn and (i'iv)am, reduces to: 

Vn = ^ ; {VNjani = ^ ; (106) 

which implies an adjoint spherical configuration with rotation parameter, 
{vN)iso = + {vN)ani = 0, as expcctcd. This is why an additional real 
rotation, due to the change Erot Ej-ot — ^E^es, is exactly balanced by 
an additional imaginary rotation, due to the change E^es E^es + ^E^es-i 
according to Eqs. (106). Then the initial shape remains unchanged. In ad- 
dition, owing to Theorem 4, the condition of stability is also expected to be 



''With the additional restriction, that streaming motions within the volume element at 

the end of major equatorial axis arc allowed provided the existence of a given flow therein 
implies a counter-flow, such that Eq. (88) is satisfied in the region of interest. 
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left unaltered. In other words, the rotating sphere (Cn = C22 < C33) is ex- 
pected to be as stable as the nonrotating sphere (Cii = C22 = Css); regardless 
from the assumed density profile, according to the results from numerical 
simulations (Meza 2002). 

4 Application: physical interpretation of the 
early Hubble sequence 

The results of the current paper allow an extension and a generalization to 
anisotropic peculiar velocity distributions, of the classical physical interpre- 
tation of the Hubble (1926) sequence (e.g., Jeans 1929, Chap. XIII, §§298- 
303). In the following, quotations are from therein. Owing to Theorem 4, 
it will suffice to restrict to solid-body rotating configurations with isotropic 
peculiar velocity distributions, but with imaginary rotation also taken into 
consideration. 

4.1 Classification of galaxies 

It is worth recalling that, before establishing their similarity to the Milky 
Way, galaxies were referred to as "great nebulae" or "nebulae" . 

Hubble finds that it is not possible to place all observed neb- 
ulae in one continuous sequence; their proper representation de- 
mands a Y-shaped diagram. (...) 

The lower half of the Y is formed by nebulae of approxi- 
mately elliptical or circular shape. These are subdivided into 
eight classes, designated £'0, El, E7, the numerical index be- 
ing the integer nearest to 10[(a — b)/a], where a and b are the 
greatest and the least diameter of the nebulae as projected on 
the sky. Thus EO consists of nearly circular nebuale (6 > 0.95a), 
while E7 consists of nebulae for which b is about 0.3a, this being 
the greatest inequality of axes observed in the "elliptical" nebu- 
lae. (...) 

The upper half of the Y-shaped diagram consists of two dis- 
tinct branches, one of which is found to contain a far larger num- 
ber of nebulae than in the other. The principal branch contains 
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the normal "spiral" nebulae, which arc characterized by a circu- 
lar nucleus from which emerge two (or occasionally more) arms of 
approximately spiral shape. (...) These nebulae are subdivided 
into three classes, designated Sa, Sb, Sc, class Sa fitting almost 
continuously on to class E7. 

The minor branch contains a special class of spirals character- 
ized by the circumstance that the spiral arms appear to emerge 
from the two ends of a straight bar-shaped or spindle-shaped 
mass. (...) 

About 97 per cent of known extra-galactic nebulae are found 
to fit into this Y-shaped classification. The remaining 3 per cent, 
are of irregular shape, and refuse to fit into the classification at 
all. (...) The irregular nebulae are distinguished by a complete 
absence of symmetry of figure and also by the absence of any 
central nucleus. (...) 

Apart from the irregular nebulae, Hubble states that, out of 
more than a thousand nebulae examined, less than a dozen re- 
fused to fit into the Y-shaped diagram at all, while in less than 
10 per cent of the cases was there any considerable doubt as to a 
proper position of a nebula in the diagram. Clearly, the Y-shaped 
diagram provides a highly satisfactory working classification. 

Thought some improvements have been made to the classification of 
galaxies (e.g., van den Bergh 1960a,b; Roberts & Haynes 1994), still we 
shall maintain Hubble (1926) diagram as a valid point of reference, owing to 
its intrinsic simplicity. 

4.2 Physical interpretation 

Obviously the proper physical interpretation of the classifica- 
tion just described is of the utmost importance to cosmogony. 

A first and most important clue is provided by the fact that 
numbers of the great nebulae are known to be in rotation. (...) 

The symmetry of figure shewn by nebulae of the E and Sa 
types is precisely such as rotation might be expected to produce, 
and this suggests an inquiry as to how far the observed figures of 
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the regular nebulae can be explained as the figures assumed by 
masses rotating under their own gravitation. (...) 

Remembering that rotation has actually been observed in a 
number of nebulae, there seem to be strong reasons for conjec- 
turing that the observed configuration of the nebulae may be ex- 
plained in general terms as the configurations of rotating masses. 

More specifically, two classes of models are considered, where density 
profiles range between the limiting cases (i) homogeneous configurations i.e. 
classical MacLaurin spheroids and Jacobi ellipsoids, and (ii) mass points sur- 
rounded by massless atmospheres i.e. Roche systems. Rigidly rotating mass 
distributions of the kind considered may attain any shape between the ex- 
treme boundaries (a) spherical i.e. nonrotating, and (b) centrifugal support 
at the end of major equatorial axis. With respect to the latter configura- 
tion, further rotation implies equatorial shedding (sequence SN) or top ma- 
jor equatorial axis streaming (sequenceSB) of matter, for axisymmetric and 
triaxial configurations, respectively. For a fixed density profile, the sequence 
of rigidly rotating configurations starts from the spherical shape and ends at 
the shape where ccntrifTigal support is first attained. 

The above interpretation suffers from two main points. First, the most 
fiattened elliptical configuration should occur when centrifugal support takes 
place at the bifurcation point, from axisymmetric to triaxial shapes. On 
the other hand, the above mentioned limiting configuration is less fiattened 
than El. Second, triaxial bodies of the kind considered cannot be figures of 
equilibrium if a great central condensation of mass is present. On the other 
hand, it is the case for barred spirals. 

Owing to Theorem 4, homeoidally striated Jacobi ellipsoids with arbitrary 
peculiar velocity field, and systematic motions reduced to rotation around a 
fixed axis, may be related to classical Jacobi ellipsoids, provided imaginary 
rotation is taken into consideration. Then the occurrence of anisotropic pe- 
culiar velocity distributions does not affect the physical interpretation under 
discussion, and the questions to be clarified are the above mentioned two. 

4.3 The E sequence within EUipsoidland 

In dealing with a physical interpretation of the early Hubble sequence, it is 
convenient to define axis ratios of intrinsic configurations in a different way. 
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Given a homcoidally striated Jacobi ellipsoid, let a, b, c, be the semiaxes, 
where a > b > c without loss of generality. Accordingly, tea ^ ^ba ^ 1; in 
addition, the minor and the major axis coincide with the rotation axis for 
oblate-like and prolate-like configurations, respectively. 

The whole range of possible configurations in the (OefeaCca) plane defines 
Ellipsoidland (the term is from Hunter & dc Zccuw 1997). EUipsoidland is 
a triangle where two sides are of unit length and an angle is right, as shown 
in Fig. 3. Oblate configurations lie on the vertical side, = 1, the spherical 
configuration is on the top, Cca = 1, and the fiat circular configuration is on 
the bottom, eca = 0. Prolate configurations lie on the inclined side, = ^ha, 
the spherical configuration is on the top, tea = Cfoa = 1, and the prolate oblong 
configuration is on the bottom, Cca = Cfoa = 0. Flat configurations lie on the 
horizontal side, — 0, the flat circular conflguration is on the right, — 1, 
and the fiat oblong configuration is on the left, = 0. It is worth noticing 
that fiat-oblong and prolate-oblong configurations are not coincident (e.g., 
Caimmi 1993; CMOS). Non fiat, non axisymmetric configurations, lie within 
Ellipsoidland. 

The axis ratio correlation, eca vs. €(,„, related to adjoint configurations 
to homeoidally striated Jacobi ellipsoids, is represented within Ellipsoidland 
in Fig. 3. The caption of symbols is the same as in Fig. 1. With regard to 
real rotation, the correlation is the known one involving classical MacLau- 
rin spheroids and Jacobi ellipsoids. With regard to imaginary rotation, the 
correlation reads eca = ^ba, owing to lack of bifurcation points. Then con- 
figurations in real and imaginary rotation branch off from the nonrotating 
spherical configuration. The former sequence proceeds down along the oblate 
side of Ellipsoidland, until the bifurcation point is attained and the curve en- 
ters Ellipsoidland, finishing when the next bifurcation point is also attained. 
The latter sequence proceeds down along the prolate side of EUipsoidland, 
until the prolate-oblong configuration is attained, never entering Ellipsoid- 
land. On the other hand, a sequence may stop earlier, when centrifugal 
support at the ends of major equatorial axis is attained, which depends on 
the density profile. In fact, it makes the sole difference between homeoidally 
striated Jacobi ellipsoids and their adjoint counterparts (for further details, 
see CMOS). 

With regard to Ellipsoidland, the observed lack of elliptical galaxies more 
fiattened than E7 translates into the inequality, eca > 0.25, the treshold 
being represented by the lower dashed horizontal line in Fig. 3. The locus of 
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Figure 3: Axis ratio correlation within Ellipsoidland, with regard to adjoint 
configurations to homeoidally striated Jacobi eUipsoids. The semiaxes are a, 
b, c, where a > b > c or ^ ^ba — 1 ) without loss of generahty. The locus of 
flat conflgurations is eca = 0, 1- The locus of oblate conflgurations 

is Cfea = 1, < < 1- The spherical conflguration is defined as Cca = ^ba = 1- 
The oblong configuration is defined as Sca = ^ba = 0. The locus of prolate 
configurations is e^a = ^ba; < e^^ < 1. Different configurations occur in the 
flat oblong limit, = 0, e^a ^ 0, and in the prolate oblong limit, — >^ 0, 
^ba ~^ 0, Eca — ^ba- A band, bounded by two horizontal dashed lines, deflnes 
class El in Hubble (1926) classification i.e. 0.25 < e^a < 0.35 for edge-on 
ellipsoids where the line of sight coincides with the direction of the minor 
equatorial axis. Other symbols as in Fig. 1. 



edge-on configurations, viewed along a direction coinciding with equatorial 
minor axis, where 0.25 < eca < 0.35, and then belonging to class E7, is 
represented in Fig. 3 as a band bounded by two dashed lines, Sca — 0.25 
and Sea = 0.35, respectively. A change in direction of the line of sight could 
project an intrinsic ellipsoid within the above mentioned band, on a class Ei, 
i < 7. 

An inspection of Fig. 3 shows that the bifurcation points from triaxial and 
axisymmetric configurations, towards pear-shaped configurations, are consis- 
tent, or marginally consistent, with the E7 band on Ellipsoidland. The above 
considerations provide a physical interpretation to the observed absence of 
elliptical galaxies (and spiral bulges) more flattened than E7, with regard 
to the oblate-like branch of the sequence. On the other hand, the above 
interpretation cannot apply to the prolate branch of the sequence, as no bi- 
furcation point has been found therein, and some kind of instability must be 
considered. 

The absence of elliptical galaxies more flattened or elongated than E7 
might be due to bending instabilities, as suggested long time ago from an- 
alytical considerations involving homogeneous (oblate and prolate) spheroids 
(Polyachenko & Shukhman 1979; Fridman & Polyachenko 1984, Vol. 1, Chap. 4, 
Sect. 3.3, see also pp. 313-322; Vol. 2, p. 159) and numerical simulations in- 
volving inhomogeneous (oblate and prolate) spheroids (Merritt & Hernquist 
1991; Merritt & Sellwood 1994). The amount of figure rotation seems to be 
unimportant to this respect (Raha et al. 1991; Merritt & Sellwood 1994). 

This conclusion is supported by recent results from A'"-body numerical 
simulations, where spherically symmetric, unstable, radially anisotropic, one- 
component 7 models were taken as starting configurations (Nipoti et al. 
2002). The related end-products, in accordance with previous results (e.g., 
Merrit & Aguilar 1985; StiaveUi & Sparke 1991), were found to be in general 
prolate systems less flattened than E7 (Nipoti et al. 2002, Fig. 1 therein). 

The theoretical explanation of the result, in terms of a dynamical bending 
instability (Merritt & Sellwood 1994), is generally recognized to explain also 
the maximum elongation of simulated nonbaryonic dark matter haloes (e.g., 
Bett et al. 2007). 
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4.4 Tidal effects from hosting dark matter haloes 

Current ACDM cosmologies, which provide a satisfactory fit to data from 
primordial nucleosynthesis and cosmic background radiation, predict galaxies 
are embedded within dark matter haloes. Then it is a natural question 
to what extent the presence of hosting dark matter haloes may affect the 
above interpretation of the early Hubble sequence. To this aim, an idealized 
situation shall be analysed. 

Let us represent elliptical galaxies and their hosting dark haloes as con- 
centric and coaxial classical Jacobi ellipsoids, one completely lying within 
the other. Accordingly, the two bodies must necessarily rotate at the same 
extent, and/or one of them (the outer in the case under discussion) has to 
be axisymmetric. It can be seen that the effect of the outer ellipsoid on the 
inner one is to shift bifurcation points towards more flattened configurations 
with respect to a massless embedding subsystem (Durisen 1978; Pacheco et 
al. 1986; Caimmi 1996a). 

More specifically, the axis ratio of the configuration at the bifurcation 
point is the solution of the transcendental equation (Caimmi 1996a): 



(A) 



5-4[(e,)3i]'^ m i-A[{e,U^A,] 



3 - 2[(e,)3i]2 ymVs 3 - 2[(e,)3i]2 (A)3 



;(107a) 



„ = § ; J,,= (^ ; ^ = ; (107b) 

Mi {ai)k yiy2y3 Pi 

where the indices, i, j, denote inner and outer ellipsoid, respectively. 

The effect of the embedding subsystem on the shape of the embedded 
configuration at the bifurcation point, is governed by the parameter: 

n= [^j]^ ; 0<fi;<+oo ; (108) 

ymys (Aja 

in the special case of massless outer ellipsoid, k, — 0, Eq. (107a) reduces to 
(65). 

Further analysis shows that the embedded configuration, at the bifurca- 
tion point, is as flattened as E7 when the parameter, k, is close to unity. 
Accordingly, the following relation holds: 

\ ; 0.25 < (q)3i < 0.35 ; (109) 
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where the shape factor ratio, {Aj)3/{Ai)s, exceeds unity if the outer ellipsoid 
is more flattened than the inner one, and vice versa. 

The generahzation of Eq. (109) to homeoidally striated Jacobi eUipsoids 
demands to restart from a generalized formulation of the virial theorem, 
where the tidal potential is also included (Brosche et al. 1983; Caimmi et al. 
1984; Caimmi & Sccco 1992). The repetition of the same procedure used in 
the current paper, yields: 

m iiy,,)uAA,h _^ . 0.25 < (e,)3i < 0.35 ; (110) 



yiy2y3 ii^i)sei (^i)3 



where {uij)tid is an additional factor in the expression of the potential tidal 
energy, related to the tidal action of the embedding subsystem on the em- 
bedded one (Caimmi 2003; CM05). 

It may safely be expected that the external boundary is less flattened than 
the internal one and more flattened than a sphere. Accordingly, the shape 
factor ratio appearing in Eqs. (109), (110), ranges as 1/2 < (Aj)3/(Aj)3 < 1 
for oblate configurations, and 1 < {Aj)3/{Ai)3 < 10/3 for prolate configura- 
tions. Then, to a first extent, (^4^)3/(^^)3 ~ i-e. the boundaries are similar 
and similarly placed ellipsoids. With this restriction, the factor, {i'ij)tid, is a 
profile factor, which may be expressed by a simple formula (Caimmi 2003): 

, "^f, n^^'-'^^ ; (111) 

where yo is the scaling radius ratio of outer to inner subsystem, and w^^^^rj) 
is an additional profile factor. 

Typical elliptical galaxies and their hosting dark haloes may safely be 
represented as homeoidally striated Jacobi ellipsoids, where the star and 
dark subsystem are described by generalized power-law density profiles of 
the kind: 

2^ /9 - 7 

^i^yo^j ; yEi^yoEj ; (112c) 
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according to Eqs. (1), and the choices 7) = (1,4,1) (Hernquist 1990) 

and (a,/9,7) = (1,3,1) (Navarro et ah 1995, 1996, 1997) are adopted to 
describe the star and the dark subsystem, respectively. The values of input 
and output parameters of the model are listed in Tab. 1, with regard to a 
ACDM cosmology where VIm = 0.3, Ha = 0.7, = O-OUbh'^, h = 2-^2. 
For further details and references, sec CM05. 

The presence of a massive halo has little influence on the location of the 
bifurcation point from axisymmetric to triaxial configurations, which is found 
to occur at an axis ratio, (esi)^^/ = 0.580088. A similar result is expected 
to occur for the bifurcation point from both axisymmetric and triaxial to 
pear-shaped configurations. Then the above mentioned points continue to 
be consistent, or marginally inconsistent, with the E7 band on EUipsoidland, 
even in presence of (typical) massive dark haloes, concerning the oblate-like 
branch of the sequence. On the contrary, it is suggested the existence of 
some kind of instability, which does not allow prolate configurations in rigid 
imaginary rotation, more elongated than E7, even in presence of (typical) 
massive dark haloes, according to recent results from iV-body simulations 
(Nipoti el al., 2002). 

4.5 Cosmological effects after decoupling 

About thirthy years ago, Thuan & Gott (1975) idealized elliptical galax- 
ies as MacLaurin spheroids, resulting from virialization after cosmological 
expansion and subsequent collapse and relaxation of their parent density 
perturbation. A generalization of the method to triaxial configurations and 
anisotropic peculiar velocity distributions, has been performed in CM05, and 
an interested reader is addressed therein for further details. Owing to Theo- 
rem 4, the case of isotropic peculiar velocity distribution, involving both real 
and imaginary rotation, can be considered without loss of generality. 

Our attention shall be limited to dark matter haloes hosting giant galax- 
ies, as massive as about lO^^m©. Accordingly, it is assumed 6rec. = 0.015 
as a typical overdensity index of the initial configuration, taken to be at 
recombination epoch, and a relaxed final configuration i.e. C = 1- 

With regard to the initial configuration, the following approximations 
hold to a good extent: (i) spherical shape; (ii) homogeneous mass distribu- 
tion; (iii) negligible (real and/or imaginary) rotation energy; (iv) negligible 
peculiar energy. The changes in mass, angular momentum, and total energy. 
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input 


value 


output 


value 




40/3 




10.38399 




10 




17.86565 


(ro)i/kpc 


3.21 


{^i)sel 


1.44444 


(ro)j/kpc 


36.20 


{^j)sel 


0.6268271 


-Rj/kpc 


42.77 


{^ij)tid 


0.3518317 


Rj /kpc 


362.04 


{l^ji)tid 


0.2217760 




8.33 


{^i)rot 


1/3 


Mj/lO^^rriQ 


91.67 


{^j)rot 


1/3 


yo 


11.29 


{Vi)rot 


1/2 


y 


8.47 


{Vj)rot 


1/2 


m 


11 


yjiest) 
yjiint) 

( ^ij ) int 
{yji)int 


-0.3955800 

-3.564028 
1.760852 
1.760852 
0.230769 



V ^sel J j 



_M 



0.531779 

0.00441624 
0.580088 



Table 1: Values of input and output parameters in modelling elliptical galax- 
ies and their hosting dark haloes as similar and similarly placed, homeoidally 
striated Jacobi ellipsoids, characterized by Hernquist (1990) and Navarro 
et al. (1995, 1996, 1997) density profiles, respectively. The parameters of 
the related ACDM cosmology have been chosen as Qm — 0.3, Qa — 0.7, 
Qb — 0.0125/i~^, h — 2~^/^. For further details on the output parameters see 
e.g., CM03. 
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respectively, during the transition from the initial to final configuration, are 
expressed by the parameters: 



/3m — 



M 



J_ 
J' 



E 
'E' 



(113) 



where the initial configuration is marked by the prime. For further details, 
see CMOS. 

Though negligible with respect to the potential and expansion energy, 
the rotation energy of the initial configuration affects the shape of the final 
configuration, as: 



1/2 



osc 



E' \ 

•pec) 



Pe\~\ -^^^rot 
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(114a) 
(114b) 

(114c) 

(114d) 

(114e) 



where different changes in mass, angular momentum, and tidal energy, are 
related to a same final shape for an assigned initial configuration, provided 
Ps does not vary. The choice, Ps = 1/3600, in particular {Pm, Pj, Pe) — 
(1,60, 1), holds for dark matter haloes hosting giant elliptical galaxies, and 
it shall be assumed here. For further details, see CM05. 

The axis ratios of the final configuration, e^i and €21, as a function of the 
parameter, K£ = S^^^/ (1 — S'^^^ — Sp^c), are plotted in Fig. 4, with regard to 
a Navarro et al. (1995,1996,1997) density profile (S = 9.20678; see CM05, 
Tab. 1 therein, for further details) in rigid rotation. Each curve is symmetric 
with regard to a vertical axis, k,£ = 0.5, where a local minimum is attained. 
The plane, (O^c^e), may be divided into three regions, namely (a) 0.163190 < 
K£ < 0.836810, where triaxial configurations occur, ranging from (e2i,e3i) = 
(1,0.582724) to (621,631) = (0.363174,0.302297), the latter related to K£ = 
0.5; (b) 0<K£< 0.163190, 0.836810 < < 1, where oblate configurations 
occur, ranging from (621,^31) = (1,0.582724) to (6215^31) = (1; 1); the latter 
related to k,£ = (bound), k,£ = 1 (unbound), respectively; (c) —00 < 
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Figure 4: Equatorial (upper curve) and meridional (lower curve) axis ratio 
of the final configuration, as a function of the parameter, K£ = S^^^/{1 — 
^'osc ~ ^pec)^ for sequences of relaxed Navarro et al. (1995, 1996, 1997) 
density profiles in rigid rotation, related to dark matter haloes hosting gi- 
ant elliptical galaxies. The curves are symmetric with respect to a vertical 
axis, K,£- = 0.5, and a minimum extremum point occurs at {k£,€2i,€^i) = 
(0.5,0.363174,0.302297), corresponding to the most flattened, oblate-like 
configuration which is allowed. Values of Kg within or outside the range, 
< Kg < 1, are related to real or imaginary rotation, respectively. Values of 
below or not below unity, are related to bound (finite extent) and unbound 
(infinite extent) configurations, respectively. Oblate-like (631 < 621 < 1) 
triaxial configurations lie between the inner, vertical dotted lines. Oblate 
(csi < £21 = 1) axisymmetric config4^tions lie between the inner and the 
corresponding outer, vertical dotted lines. Prolate (e2i = 1 < 631) axisym- 
metric configurations lie outside the outer, vertical dotted lines. The equa- 
torial axis ratio is e = 621. The polar axis ratio is e = 631 for oblate and 
e = 613 for prolate configurations. The horizontal dashed lines define class 
E7 in Hubble (1926) classification i.e. 0.25 < e < 0.35 for edge-on ellipsoids 
where the line of sight coincides with the direction of the minor equatorial 




K-£ < 0, 1 < Kg < +00, where prolate configurations occur, ranging from 
(e2i,e3i) = (1,1) to (e2i,e3i) = (1,+cxd), tlie latter related to k.£ — >• —00 
(bound), K£ — > +00 (unbound), respectively. It is worth recalling that oblate- 
like (esi < 621 < 1) and prolate-like (621 < 1 < £31) configurations, are related 
to real and imaginary rotation, respectively. 

An inspection of Fig. 4 shows that, in the case under discussion, cosmolog- 
ical effects due to expansion prevent dark matter haloes in real rotation from 
being more flattened than E7 (dashed horizontal band) . The same holds for 
embedded giant elliptical galaxies, provided the related shapes may safely be 
thought of as similar and similarly placed. On the other hand, an arbitrary 
flattening can be attained by dark matter haloes in imaginary rotation, un- 
less some kind of instability occurs, which similarly prevents configurations 
more elongated than E7. 

As outlined in Subsect. 4.3, bending instabilities have been suggested long 
time ago as a viable mechanism to this respect (Polyachenko & Shukhman 
1979; Fridman & Polyachenko 1984, Vol. 1, Chap. 4, Sect. 3.3, see also pp. 313- 
322; Vol.2, p. 159; Merritt & Hernquist 1991; Merritt & Sellwood 1994). 

The above results depend on both density profile and rotation velocity 
profile of the final configuration, via the coefficient, c, appearing in Eq. (114a). 
Owing to Eqs. (7) and (17c), Eq. (114c) translates into: 



'ram 



2C-1 C-ic 
2C C 



:il5a) 



(^s' = fis^^^^ ; (115b) 

where different changes in mass, angular momentum, total energy, density 
profile, and rotation velocity profile, are related to the same final shape for 
an assigned final configuration, provided f3£i does not vary. The choice, jSg' = 
6.007519-10-^ in particular (/?m, Pj, Pe, ^seh ^ram) = (1, 60, 1, 0.610045, 17.20783), 
holds for Fig. 4. 



5 Conclusion 

With regard to homeoidally striated Jacobi ellipsoids (CM05), a unified the- 
ory of systematically rotating and peculiar motions have been developed. 



45 



where both real and imaginary rotation (around a fixed axis) have been con- 
sidered. The effect of residual motion excess along the equatorial plane or the 
rotation axis, has been shown to be equivalent to an additional, real or imag- 
inary rotation, respectively. Then it has been realized that a homeoidally 
striated Jacobi ellipsoid with assigned velocity distribution, always admits 
an adjoint configuration i.e. a classical Jacobi ellipsoid of equal mass and 
axes. In addition, further constraints have been established on the amount 
of residual velocity anisotropy along the principal axes, for triaxial configu- 
rations. Special effort has been devoted to investigating sequences of virial 
equilibrium configurations in terms of normalized parameters, including the 
effects of both density and velocity profile. In particular, it has been shown 
that bifurcation points from axisymmetric to triaxial configurations occur as 
in classical Jacobi eUipsoids, contrary to earher results (Wiegandt 1982a,b; 
CM05). The reasons of the above mentioned discrepancy have also been 
explained. The occurrence of centrifugal support at the ends of equato- 
rial major semiaxis, has briefiy been discussed. An interpretation of recent 
results from numerical simulations on stability (Meza 2002), has been pro- 
vided in the light of the model. A physical interpretation of the early Hubble 
sequence has shortly been reviewed and discussed from the standpoint of 
the model. More specifically, (i) elliptical galaxies have been considered as 
isolated systems, and an allowed region within Ellipsoidland (Hunter & de 
Zeeuw 1997), related to the occurrence of bifurcation points from ellipsoidal 
to pear-shaped configurations, has been shown to be consistent with ob- 
servations; (ii) elliptical galaxies have been considered as embedded within 
dark matter haloes and, under reasonable assumptions, it has been shown 
that tidal effects from embedding haloes have little infiuence on the above 
mentioned results; (iii) dark matter haloes and hosted elliptical galaxies, ide- 
alized as a single homeoidally striated Jacobi ellipsoid, have been considered 
in connection with the cosmological transition from expansion to relaxation, 
by generalizing an earlier model (Thuan & Gott 1975), and the existence 
of a lower limit to the fiattening of relaxed (oblate-like) configurations, has 
been established. On the other hand, no lower limit has been found to the 
elongation of relaxed (prolate-like) configurations, and the observed lack of 
elliptical galaxies more elongated than E7 has needed a different physical 
interpretation such as the occurrence of bending instabilities (Polyachenko 
& Shukhman 1979; Merritt & Hernquist 1991). 
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A Some properties of ellipsoid shape factors 

Ellipsoid shape factors, Ai, A2, A^, obey the following inequalities (Caimmi 
1996a): 

a^Ai > a'^A2 > a'^A^ ; Oi > 0,2 > ; n > 2 ; (116a) 
<Ai < a!^A2 < a'^A^ ; Oi > as > 03 ; n < 1 ; (116b) 

where inequalities (116a) and (116b), the latter restricted to n < 0, come 
from analytical considerations involving the explicit expression of the shape 
factors (e.g., MacMillan 1930, Chap. 11, §33; Ai = a'^, A2 = A3 = 7-2, 
therein), and inequality (116b), restricted to < n < 1, results from an 
obvious generalization of a proof by Pacheco et al. (1989). 
Using Eqs. (6) and (7) yields: 

#^ = 3^; ^ = 1.2; (117) 

and, owing to inequality (116a): 

Ag>e%A3 ; a, > 03 ; ? = 1,2 ; (118) 

which implies Sgg > S33 for oblate-like configurations (a^ > 03), and Sgg < 
Sss for prolate-like configurations (og < 03). 

In the limit of axisymmetric configurations, €21 = 1, £31 = e, Ai = A2 = a, 
A3 — 7, and the following relations hold (e.g., Caimmi 1991, 1993): 

CK TT 

lima^O ; lim7 = 2 ; lim - = - ; (119) 

e— >0 e— >0 e— >0 e 2 

lim a — 1 ; lim 7 = lim (57) = ; lim (e^7) = -|-oo ; (120) 
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lim 



7 



a 



*i 1 
da 

^ ~~ 
d7 

V 



2 
5 



1 - e2 
1 



:i + 2e^ 



-2e 



1 -e^ 



(l + 2e2)^-^ 
a — e^7 ; 



(121) 
(122) 

(123) 

(124) 



and keeping in mind the general property (e.g., Chandrasekhar 1969, Chap. 3, 
§17): 

^1 + ^2 + ^3 = 2 ; (125) 

it can be seen that the first derivative of the rotation parameter, is null 
provided the transcendental equation: 

(126) 



a 



l + 8e2 



is satisfied. One solution is found to exist, which is related to oblate configu- 
rations. It corresponds to an extremum point where the function attains its 
maximum value (e.g., Chandrasekhar 1969, Chap. 5, §32), as represented in 
Fig.l. 



B Wiegandt criterion for bifurcation 

Given a coUisionless, self-gravitating fluid in rigid rotation, where no internal 
energy transport occurs and the residual velocity is constant on the boundary, 
an upper limit for the point of bifurcation is (Wiegandt 1982a,b): 

^^'/ll = l"l2;12 ; (127) 

where Vpq.^rs is the "super-matrix" : 

f dVrs 
Vpqyrs^ / p{xi,X'2.,Xz)Xp — S ] (128) 

J S 

V..(x„ .2, xs) ^Gj^ j^^^ _ ^ _ ^ _ ^,^3/, d S ei29) 

which is expressed in terms of a generalized potential, Vrsi the integra- 
tions are carried over the whole volume, 5, of the system. 
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In the special case of homeoidally striated Jacobi ellipsoids, Eq. (127) 
identifies the exact point of bifurcation (Wiegandt 1982a,b) and the following 
relations hold (Wiegandt 1982b): 

%>g;pg — 1 {Esei)pp (130) 



alA,,^-^ ; (131) 
pg 

where the products, a^Ap^, are shape factors which depend on the axis ratios 
only, similarly to Ap, and symmetry with respect to the indices holds, Apq = 
Agp. In addition, Eq. (131) has been deduced from Chandrasekhar (1969, 
Chap. 3, §21), Eq. (107) therein. 

Though Wiegandt (1982b) analysis is restricted to binomial density pro- 
files (e.g., Perek 1962; Chandrasekhar 1969, Chap. 3, §20; Caimmi 1993), 
still it may be generalized to any kind of cored density profiles, defined as: 

PwiO-Pcmw) ; /(0) = 1 ; Pc = p(0) ; (132a) 
^w = ^ ; o<Cw<i ; ^w^i ; (i32b) 

where the scaling radius and the scaling density are chosen to be equal to the 
radius, R, and the central density, pc, respectively. The following relations: 

P = —Pw ; C = '^Cw ; (133) 

Pc 

allow conversion from Eqs. (1) to Eqs. (132) and vice versa. 

Using Eqs. (132), the mass, the inertia tensor, and the potential-energy 
tensor, take the equivalent expression: 

M^{ur„as)wMc ; (134a) 
Mc = YPca?e2ie3i ; (134b) 



Ipq — 5pqepieqia\M[i'inr)w ] (135) 

{Esel)pq — —{^sel)w {Bsel)pq ] (136) 
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and the comparison with Eqs. (2), (3), (4), (8), yields: 

{l^mas)w — — ™i^mas ; {'^inr)w — ^inr '■, {'^sel)w — ^sel ] (137) 

finally, using Eqs. (6) and (134)-(137), the potential-energy tensor takes the 
equivalent form: 

{Esei)pq = -kTrGpJpqAp ; (138a) 

^^37^^ ' ^^^^^^ 

being k, by definition, a profile factor. 

Owing to Eqs. (8), (18), (20), (34), (35), (138b), the normalized rotation 
parameter, defined by Eq. (48), may be expressed as: 

which coincides with the parameter, 0, defined in Wiegandt (1982b), as 
shown in Caimmi (1996b). Accordingly, Eq. (49) coincides with its Wiegandt 
(1982b) counterpart, Eq. (52) therein. 

The combination of Eqs. (127), (130), (138), (139), yields: 

VN^Ai-ajAu ; (140) 
and the comparison with Eq. (49) reads: 

alAu = ^el,As ; (141) 

^,33 

in the limit of axisymmetric configurations, 02 ~^ cli, ^12 ~^ ^11, the left- 
hand side of Eq. (141) takes the expression (e.g., Caimmi 1993): 

aF-^„,«2^- = 1(3^,-44^) ; (142) 



and the combination of Eqs. (141), (142), yields: 



-31^3 

A, 



5-46^, ^ 4(l-e^JCn-C33 
3 — 2631 3 — 2e3i (^33 



£21 = 1 ; (143) 
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which is an exphcit expression of Wiegandt (1982a,b) criterion for bifurca- 
tion, with regard to homeoidally striated Jacobi elhpsoids in rigid rotation. 
In the hmit of isotropical residual velocity distribution, di = C22 = Csa? 
Eq. (143) reduces to (65), as expected. 

The condition, defined by Eq. (143), has to be compared with its counter- 
part deduced in the current paper, under the same assumption i.e. (n = (22 
also for triaxial configurations. The latter may be deduced from Eq. (84), 
using the relation (e.g., Caimmi 1966a): 

2 A2-A l/„. 2 As-AA ..... 

-"'"T^ = 4(,'^^-^-^^j ; (144) 

and the result is again Eq. (143). Then Wiegandt (1982b) criterion for bifur- 
cation can be deduced from the current theory, provided isotropic residual 
velocity distribution is assumed along the equatorial plane, also for triaxial 
configurations. On the other hand, the above assumption is in contradiction 
with Eqs. (63) and then it cannot be accepted. A criterion for bifurcation, 
consistent with Eqs. (63), is expressed by Eq. (65). 

The origin of the discrepancy could be the following. Wiegandt (1982a) 
assumption (ii) implies Eq. (12) therein which, after integration, yields isotropic 
residual velocity distribution along the equatorial plane i.e. Cii = (22- On 
the other hand, an isotropic residual velocity distribution along the equato- 
rial plane is related to a zero-th order approximation in some special series 
developments [Marochnik 1967; Eq. (19) therein]. The additional require- 
ment of axial symmetry makes the above mentioned relation be exact, but 
in presence of triplanar symmetry it has to be considered as a zero-th order 
approximation. 

In the light of the results discussed in the current subsection, the anal- 
ogy between the behaviour of collisional and collisionless self-gravitating fiu- 
ids, subjected to the restrictions (Wiegandt 1982a,b): (i) rigid rotation; (ii) 
constant residual velocity on the boundary; (iii) absence of internal energy 
transports; appears to be complete. 
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